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IGNITION  OF  FLAMELETS  BEHIND  INCIDENT  SHOCK  WAVES 
AND  THE  TRANSITION  TO  DETONATION 

I.  INTRODUCTION 

Shock-induced  ignition  in  gaseous  hydrogen-air  mixtures  may  occur  in  one 
of  two  distinct  modes,  depending  on  the  thermodynamic  state  in  the  shocked 
material.  At  lower  temperatures  the  ignition  is  weak,  or  mild,  with  the 
gradual  development  of  the  gas  dynamic  explosion.  At  higher  tenperatures  it 
is  strong,  or  sharp,  with  a.%  abrupt  appearance  of  a  secondary  shock  induced 
by  the  explosive  reaction.  Furthermore  at  low  tenperatures  the  formation  of 
distinct  flame  kernels  appears  to  be  an  essential  precursor  of  ignition. 

Meyer  and  Cppenheim  (1971a)  point  out  that  the  intrinsically  turbulent  flow 
field  behind  a  reflected  shock  in  a  shock  tube  results  in  a  nonuniform 
teuperature  distribution  which  creates  distinct  reaction  centers  and  leads  to 
weak  ignition.  Although  the  nonuniformity  of  temperature  caused  by 
turbulence  is  one  mechanism  which  gives  rise  to  reaction  centers,  the 
nonsteadiness  in  the  velocity  of  the  causal  shock  can  also  produce  them. 

This  latter  effect  has  been  convincingly  demonstrated  by  Strehlow  et  al. 
(1967)  in  their  studies  of  shock  propagation  in  a  slowly  converging  channel. 
It  has  also  been  shown  that  when  an  incident  shock  in  a  uniform  shock  tube 
begins  to  accelerate,  hot  spots  occur  before  the  transition  to  detonation 
(Bazhenova  and  Soloukhin,  1959;  Edwards  et  al. ,  1981). 

In  this  paper  the  reactive  flow  behind  an  incident  shock  wave  in  a 
hydrogen-air  mixture  is  sinulated  using  a  one -dimensional,  time -dependent 
numerical  model  which  combines  a  description  of  the  fluid  dynamics  and 
detailed  chemical  kinetics.  In  the  simulations,  the  pressure  ratio  across 

the  diaphragm  and  the  driven  gas  mixture  are  chosen  so  that  the  thermodynamic 
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state  behind  the  incident  shock  is  in  the  veak  ignition  regime  (Oran  and 
Boris,  1981b).  It  is  shown  that  a  small  amount  of  energy  released  in  the 
shocked  gas  (vhich  might  occur  due  to  density,  temperature  or  stoichiometric 
fluctuations)  can  be  the  origin  of  pressure  waves  which  accelerate  the  shock 
front.  Once  the  shock  is  accelerated,  the  tesqperature  of  the  gases  now 
passing  through  the  shock  is  higher  than  the  raised  temperature  created  by 
the  original  shock.  In  the  veak  ignition  regime,  such  an  increased 
temperature  can  result  in  a  significant  reduction  in  the  induction  time 
(Meyer  and  Oppenheim.  1971b;  Oran  and  Boris,  198Q.b).  The  simulations 
presented  here  show  how  this  leads  to  the  formation  of  reactive  centers  in 
the  newly  shocked  material  where  reaction  progresses  at  a  more  rapid  rate 
than  in  the  previously  shocked  material.  The  formation  of  a  hot  spot  due  to 
energy  release  at  one  of  the  reactive  centers  and  the  subsequent  development 
of  a  pair  of  flame  lets  or  reaction  waves  from  the  hot  spot  are  studied  using 
the  numerical  simulations.  The  results  of  the  simulations  have  also  been 
compared  to  experimental  observations  (Bazhenova  and  Soloukhin,  1959;  ttrtlev 
and  Oppenheim,  1966,  1967;  Strehlov  et  al. ,  I96T;  Edwards  et  al. ,  1981). 

II.  THE  NUMERICAL  MODEL 

The  one-dimensional  reactive  shock  model  (Oran  et  al.,  1979;  Oran  and 
Boris,  1981a)  used  to  perform  the  calculations  described  below  solves  the 
time-dependent  conservation  equations  (Williams,  1965,  Oran  and  Boris,  1981a) 
for  mass,  momentum  and  energy  coupled  to  the  equations  describing  the 
chemical  kinetics.  The  model  uses  an  explicit,  Eulerian  finite  difference 
formulation  with  a  sliding  rezone  capability  to  provide  resolution  around 
moving  gradients.  The  solutions  of  the  equations  describing  the  fluid 
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dynamics  and  the  chemistry  of  the  problem  are  coupled  using  time-step 
splitting  techniques  (Oran  and  Boris,  1981a). 

The  convective  transport  terms  in  the  conservation  equations  are  solved 
using  one  variant  of  the  Plux-Corrected  Transport  (FCT)  method  (Boris  and 
Book,  1976;  Boris,  1976).  This  is  a  conservative,  raonotonic  algorithm  with 
fourth-order  phase  accuracy  and  does  not  require  artificial  viscosity  to 
stabilize  shocks.  The  Flux  Correction  procedure  itself  ensures  that  the 
shocks  are  one  or  two  zones  vide  and  have  maximal  resolution.  The  ordinary 
differential  equations  describing  the  chemical  kinetics  are  solved  using 
VSAIM,  a  vectorized  version  of  the  selected  asynptotic  integration  method 
employed  in  CHIMEQ  (Young  and  Boris,  1977;  Young,  1980).  This  algorithm 
identifies  the  stiff  equations  for  treatment  with  a  stiffly  stable  method. 

The  remaining  equations  are  solved  with  a  standard  classical  method.  The 
algorithm  has  been  specially  optimized  for  use  in  conjunction  with  fluid 
dynamic  models. 

The  chemical  kinetics  rate  scheme  used  is  given  in  Table  X.  It  consists 
of  about  fifty  rates  relating  the  species  H^,  02,  H,  0,  OH,  HQ2,  H20  and  HjO^ 
and  has  been  extensively  tested  against  experimental  data  (Oran  et  al. ,  198l» 
Burks  and  Oran,  1980).  Burks  and  Oran  (1980)  shoved  that  the  results 
computed  with  the  scheme  compared  very  veil  vith  experimentally  observed 
induction  times,  second  explosion  limits  and  the  temporal  behavior  of 
reactive  species.  Oran  et  al.  (1981)  have  shovn  that  the  scheme  gives  good 
results  vhen  coupled  vith  a  fluid  dynamic  model  in  the  simulation  of  the 
conditions  behind  a  reflected  shock.  Heats  of  formation  and  enthalpies  have 
been  taken  from  the  JA2IAF  tables  (I97l)» 


3 


For  the  calculations  performed  in  this  paper,  the  timescales  under 
consideration  are  short  and  therefore  the  diffusive  transport  processes, 
thermal  conduction  and  molecular  diffusion,  have  negligible  effect.  The 
effects  of  these  processes  have  not  been  considered  although  they  are  part  of 
the  general  numerical  model.  Although  the  geometry  may  be  either  cartesian, 
cylindrical,  spherical  or  some  generalized  co-ordinate,  the  simulations 
presented  below  are  in  cartesian  geometry. 

Hie  detailed  simulations  discussed  in  this  paper  require  that  we  model 
relatively  long  systems  (on  the  order  of  meters)  while  we  simultaneously 
maintain  high  accuracy  around  steep  gradients  such  as  the  shock  front  and  the 
contact  surface.  The  rather  sophisticated  adaptive  gridding  method  developed 
for  this  purpose  is  shown  schematically  in  Figure  1.  There  are  two  finely 
gridded  regions:  One  surrounding  the  shock  wave  and  the  other  surrounding 
the  contact  surface.  The  fine-zoned  region  around  the  contact  surface  moves 
with  the  contact  surface  at  the  fluid  velocity,  and  so  this  part  of  the 
calculations  is  essentially  Lagrangian.  The  region  around  the  shock  front 
moves  with  the  front.  Bach  of  these  finely  gridded  regions  may  have  a 
different  minimum  computational  cell  size.  The  computational  cells  in  the 
regions  ahead  of  the  shock  wave  and  behind  the  contact  surface  change 
exponentially  in  size  from  the  smallest  near  the  shock  wave  or  the  contact 
surface  to  the  largest  at  the  walls.  Care  is  taken  that  the  transition  in 
the  cell  sizes  is  smooth.  For  the  results  presented  in  this  paper  a  total  of 
200  computational  cells  are  used  to  describe  the  shock  tube  and  the  cell 
sizes  varied  from  0.1  cm  around  the  shock  to  over  50  cm  near  the  shock  tube 
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III.  WEAK  IGNITION  BEHIND  AN  INCIDENT  SHOCK 

The  numerical  model  described  in  the  preceding  section  was  used  to  study 
weak  ignition  behind  planar  incident  shock  waves.  The  system  parameters,  the 
initial  conditions  and  the  tenperature  and  pressure  behind  the  incident  shock 
have  all  been  summarized  in  Table  II.  The  chemical  induction  time,  for  the 
conditions  described  by  the  incident  shock,  is  about  2000  "s.  However,  since 
the  thermodynamic  state  behind  the  shock  is  in  the  weak  ignition  regime, 
ignition  may  occur  raich  earlier  than  2000  *'s  due  to  tenperature,  pressure,  or 
density  fluctuations. 

Figure  2  is  a  position-time  diagram  of  the  events  occurring  in  the  shock 
tube  simulation.  The  trajectory  of  the  shock  front  is  labelled  S  and  that  of 
the  contact  surface  is  labelled  CS.  Except  for  small  variations  (which  are 
examined  in  detail  below)  the  shock  travels  at  a  nearly  constant  velocity, 

1.4  x  io 5  cm/s,  until  the  reaction  wave  formed  between  the  contact  surface 
and  the  shock  front  reaches  it.  At  this  time  the  velocity  rises  quickly  to 
3.24  x  io5  cm/s.  It  then  gradually  decreases  towards  the  Chapman-Jouguet 
detonation  velocity.  Five  different  regimes  have  also  been  identified  on  the 
diagram.  They  are  (a)  pre-ignition  regime,  (b)  quasi-steady  shock-reaction 
complex,  (c)  formation  of  reactive  centers,  (d)  hot  spot  formation  leading  to 
an  overdriven  detonation  and  (e)  detonation  relaxation.  Similar  regions  have 
been  identified  by  Edwards  et  al.  (1981)  in  their  shock  tube  experiments. 
These  regions  are  examined  in  detail  below. 

As  the  shock  travels  at  a  nearly  constant  velocity  into  the  hydrogen-air 
mixture,  the  temperature  and  pressure  of  the  mixture  are  raised  to  a  near¬ 
constant  value.  Reactions  in  the  shock  heated  gas  first  occur  near  the 
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contact  surface  since  the  temperature  has  been  high  for  the  longest  time 
here,  as  has  been  observed  by  Urtiew  and  Oppenheim  (196?).  In  Figure  3,  the 
spatial  variation  of  the  temperature  and  the  OH  mole  fraction  between  the 
shock  and  the  contact  surface  are  shown  at  135  us  after  the  bursting  of  the 
diaphragm.  The  OH  mole  fraction  attains  a  maximum  value  near  the  contact 
surface  and  decreases  across  the  system  towards  the  shock  front.  This  is 
because  the  gas  mixture  at  different  locations  has  been  at  the  higher 
temperature  for  different  durations. 

Formation  of  Reactive  Centers  and  Hot  Spots 

Small  pressure  disturbances  occur  due  to  the  energy  released  in  the 
reactions  near  the  contact  surface.  These  pressure  disturbances  travel 
forward  at  a  velocity  (in  the  laboratory  frame  of  reference)  of  1.8  x  105 
cm/s,  which  is  the  sum  of  the  sonic  and  particle  velocities  behind  the  shock. 
When  they  reach  the  shock  they  accelerate  the  shock  slightly  resulting  in  a 
temperature  increase  behind  the  shock.  This  can  be  observed  in  Figure  U 
where  the  spatial  distribution  of  the  teqperature  between  the  shock  and  the 
contact  surface  is  shown  at  four  different  times.  The  nonuniform  temperature 
distributions  can  be  explained  by  the  following  sequence  of  events:  pressure 
disturbances  originating  at  different  times  near  the  contact  surface  reach 
the  incident  shock  at  different  times  and  each  successive  pressure 
disturbance  meets  a  shock  of  slightly  different  strength  because  the  shock 
has  already  been  accelerated  slightly  due  to  previous  pressure  disturbances. 
We  know  that  in  the  weak  ignition  regime,  the  induction  time  is  very 
sensitive  to  perturbations  in  temperature  and  pressure  (Oran  and  Boris, 
1981b).  This  sensitivity,  in  fact,  is  what  produces  the  two  peaks  in  OH  mole 


fraction  distribution  (Figure  5A)  at  3^8  us  after  the  bursting  of  the 
diaphragm.  The  first  peak  is  due  to  the  development  of  reactions  near  the 
contact  surface  and  occurs  in  the  same  reactive  center  (see  Figure  2)  which 
was  observed  earlier.  How  a  second  reactive  center  has  formed  closer  to  the 
shock  front  due  to  the  higher  temperature  created  by  the  acceleration  of  the 
shock.  Again  we  observe  that  the  reactive  center  occurs  in  that  part  of  the 
shock  heated  gas  which  was  at  the  elevated  temperature  for  the  longest  period 
of  time.  Ahead  of  the  seconu  reactive  center  the  rapid  decrease  in  the 
radical  mole  fraction  is  halted  by  a  further  increase  in  the  temperature 
behind  the  shock  and  a  third  reactive  center  is  forming.  By  UU5  us  (Figure 
5B),  a  third  reactive  center  has  developed  and  the  radical  mole  fraction  is 
rising  most  rapidly  at  this  center.  The  temperature  around  this  reactive 
center  is  significantly  (about  Uo°  K)  greater  than  that  behind  the  initial 
incident  shock  wave.  In  the  weak  ignition  regime,  such  changes  in 
temperature  result  in  a  substantial  reduction  in  the  induction  time. 

Therefore  reaction  progresses  at  this  reactive  center  at  a  much  more  rapid 
rate  than  at  any  other  location  in  the  system.  This  results  in  a  "hot  spot" 
in  the  system.  The  development  of  this  hot  spot  can  be  seen  in  Figure  6 
where  the  temperature  distribution  between  the  contact  surface  and  the  shock 
front  has  been  shown  for  four  different  times.  The  tenperature  increases  by 
more  than  200°  in  30  us. 

Transition  to  Detonation 

The  events  occurring  after  the  hot  spot  formation  can  be  seen  in  Figure 
7  where  a  segment  of  the  position-time  diagram  (Figure  2)  is  presented  in 
greater  detail.  The  hot  spot  (marked  "A"  in  the  figure)  travels  with  the 
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fluid  and  continuously  releases  energy  into  it.  This  results  in  "flamelets" 
or  "reaction  waves"  ("B")  which  travel  with  respect  to  the  fluid.  Qie 
reaction  wave  moves  forward  with  the  fluid  while  another  moves  against  the 
fluid  towards  the  driver  section.  These  reaction  waves  initially  move 
subsonically  with  respect  to  the  fluid  but  are  soon  accelerated  into  steep 
detonation  type  waves  which  move  supersonically  with  respect  to  the  fluid. 
When  the  forward  moving  reaction  wave  reaches  the  shock  front,  the  shock 
velocity  abruptly  increases  to  a  high  value.  The  shock-reaction  complex  then 
moves  as  a  strong  overdriven  detonation  wave  which  decelerates  towards  the 
Chapman-Jouguet  detonation  velocity. 

The  temperature,  pressure  and  velocity  distributions  across  the  system 
at  a  particular  time  after  the  hot  spot  has  formed  are  shown  in  Figure  8a. 

By  this  time  the  energy  release  has  resulted  in  a  noticeable  pressure  rise. 
This  occurs  because  the  energy  release  is  occurring  at  nearly  constant  volume 
conditions.  The  energy  release  at  the  hot  spot  causes  a  series  of  minute 
pressure  pulses  to  propagate  both  forward  and  backward,  each  a  little 
stronger  than  the  previous  one.  A  series  of  pressure  pulses  are  produced 
since  the  energy  release  occurs  over  a  period  of  time  and  is  determined  at 
each  time  from  the  detailed  chemical  kinetic  interactions  among  the  various 
species.  These  pressure  pulses  coalesce  to  form  steepening  pressure  waves 
propagating  into  the  shocked  mixture  as  seen  in  Figure  8b.  The  time  history 
of  these  reaction  waves  has  been  shown  in  a  series  of  figures  (Fig.  8b  - 
Fig.  10b)  which  cover  the  time  period  between  6 09  us  and  688  us. 

Let  us  first  look  at  the  development  of  the  reaction  wave  which  moves 
forward  into  the  shock-heated  gas  mixture.  The  velocity  of  the  forward 


8 


moving  wave,  the  fluid  velocities  on  both  sides  of  the  wave  and  the  speed  of 
sound  in  the  gas  mixtures  on  both  sides  of  the  wave  are  given  in  liable  III 
for  a  series  of  times.  At  609  Vis,  the  forward  moving  reaction  wave  is 
supersonic.  The  velocity  of  the  fluid  with  respect  to  the  wave  decreases 
across  the  wave  but  is  supersonic  on  either  side.  Furthermore  the  pressure 
rise  across  the  wave  is  just  over  a  factor  of  two,  which  is  moderate.  Thus 
the  reaction  wave  behaves  like  a  weak  detonation  wave  at  this  time.  later, 
as  seen  in  Figure  9a,  the  pressure  rise  across  the  wave  has  increased  and  the 
wave  is  also  travelling  faster.  The  velocity  of  the  fluid  with  respect  to 
the  wave  is  still  supersonic  on  both  sides  of  the  wave.  The  pressure  rise 
across  the  wave  continues  to  increase  and  at  the  time  corresponding  to  Figure 
9b,  the  fluid  velocity  behind  the  wave  is  nearly  sonic.  The  weak  detonation 
seems  to  be  transitioning  into  a  strong  detonation.  It  does  so  later  (Figure 
9c)  when  the  fluid  velocity  changes  from  supersonic  to  subsonic  across  the 
reaction  wave.  The  pressure  rise  across  the  wave  is  also  larger  now.  The 
observed  acceleration  of  the  forward  moving  wave  into  a  strong  detonation  is 
due  to  the  nonlinear  interaction  between  chemical  kinetics  and  fluid 
dynamics.  When  the  forward  moving  wave  moves  into  the  previously  shocked 
material  there  is  a  large  pressure  and  temperature  rise  across  it  since  it  is 
a  strong  compressive  wave.  This  increase  in  the  pressure  and  the  temperature 
reduces  the  induction  time  of  the  material  which  crossed  the  wave.  Energy 
release  in  this  newly  re-conpressed  material  accelerates  the  forward  moving 
wave  further,  and  this  cycle  is  repeated  until  the  forward  moving  wave 
reaches  the  incident  shock  wave  (Figure  10a).  By  this  time  the  pressure 
spike  behind  the  reaction  wave  has  risen  to  6.6  x  io6  dynes/cm2.  Because  the 


9 


reaction  wave  is  moving  towards  the  incident  shock  at  a  speed  greater  them 
the  local  speed  of  sound  in  the  shock-heated  gas,  there  is  no  advance  warning 
of  the  over-pressure  region  until  the  main  spike  physically  arrives  at  the 
shock  front.  When  the  reaction  wave  does  coalesce  with  the  shock,  there  is  a 
very  rapid  increase  in  the  shock  speed  (Figure  7)  and  the  shock -react ion  wave 
complex  moves  as  an  overdriven  strong  detonation  wave.  Due  to  the  high  over¬ 
pressure  that  was  associated  with  the  reaction  front,  the  detonation  wave 
overshoots  the  Chapman-Jouguet  value  by  a  substantial  amount.  However,  the 
overshoot  cannot  be  sustained  by  the  reaction  and  the  overdriven  detonation 
gradually  relaxes  towards  a  Chapman-Jouguet  wave.  By  688  ms  (Figure  10b)  the 
pressure  spike  behind  the  detonation  wave  has  decreased  to  4.8  *  io6 
dynes/cm2  and  the  wave  is  moving  faster  than  the  C-J  velocity  by  only  3 6%. 

In  Figure  10b  we  also  observe  a  small  amplitude  pressure  wave  moving  into  the 
detonation  products.  This  pressure  wave  was  formed  when  the  reaction  wave 
interacted  with  the  shock  front. 

Let  us  now  look  at  the  reaction  wave  which  moves  backward  towards  the 
driver  section.  From  Figure  7  and  Thble  IV  we  see  that  this  wave  moves  at 
nearly  the  sonic  speed.  The  fluid  velocity  with  respect  to  the  wave  i3 
supersonic  on  either  side  of  the  wave.  The  pressure  rise  and  the 
acceleration  of  the  wave  are  slower  since  the  wave  is  moving  against  the 
fluid.  The  pressure  increase  behind  the  wave  is  broader  (Figure  9c)  and 
smaller  than  that  of  the  forward  moving  wave.  It  continues  to  propagate  like 
a  weak  detonation  wave  until  it  interacts  with  the  contact  surface  (Figure 
10a).  This  interaction  produces  a  pressure  pulse  which  travels  into  the 
helium  driver  gas.  In  Figure  10b  we  see  that  this  pressure  pulse  has 
produced  a  slight  temperature  increase  in  the  helium. 


IV.  DISCUSSION  AND  CONCLUSIONS 


In  this  paper  the  reactive  flow  behind  an  incident  shock  wave  has  been 
studied  using  a  time -dependent  numerical  model  which  includes  both  detailed 
chemical  kinetics  and  one-dimensional  fluid  dynamics.  The  numerical 
simulations  show  that  the  incident  shock  initially  travels  at  a  steady  speed 
leaving  behind  it  material  in  what  has  been  called  a  "quasi -steady  reaction 
complex".  From  their  extensive  experimental  studies  on  detonation,  Lee  et 
al.  (1976,  1977)  have  concluded  that  such  regimes  are  universal  predetonation 
phenomena.  This  regime  has  also  been  observed  in  the  incident  shock  tube 
experiments  of  Edwards  et  al.  (1981). 

The  calculations  have  also  shown  the  acceleration  of  the  shock  due  to 
pressure  waves  created  by  energy  release  in  the  shock  heated  gas  mixture,  a 
phenomenon  also  observed  by  Edwards  et  al.  (19SL).  Shock  acceleration  raises 
the  temperature  of  the  gases  passing  through  the  shock  and  in  the  weak 
ignition  regime,  this  results  in  a  significant  reduction  of  the  induction 
time  (Meyer  and  Oppenheira,  1971b;  Oran  and  Boris,  1981b).  This  leads  to  the 
formation  of  reactive  centers  where  reaction  progresses  at  a  more  rapid  rate 
than  in  the  previously  shocked  material.  The  development  of  a  hot  spot  due 
to  energy  release  at  one  of  the  reactive  centers  has  been  shown  in  Figure  6. 
The  presence  of  such  hot  spots  has  been  observed  earlier  in  shock  tube 
experiments  (Strehlow  et  al. ,  1967,  Edwards  et  al. ,  1981). 

Ehergy  release  at  a  hot  spot  causes  a  pair  of  flaroelets  or  reaction 
waves,  one  propagating  ahead  into  the  shock  heated  gas  mixture  and  the  other 
propagating  back  towards  the  contact  surface.  These  flamelets  initially 
propagate  at  a  subsonic  speed  with  respect  to  the  fluid.  Such  flamelets  have 
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been  observed  by  Bazhenova  and  Soloukhin  (1939)  in  their  incident  shock  tube 
experiments.  Energy  release  behind  these  flamelets  causes  pressure  waves 
which  accelerate  the  flamelets  into  detonation  type  waves.  The  reaction  wave 
moving  backward  against  the  fluid  accelerates  slowly  and  travels  at  nearly 
the  sonic  speed  till  it  reaches  the  contact  surface.  However  the  forward 
moving  reaction  wave  transitions  into  a  strong  detonation  wave  even  before  it 
reaches  the  incident  shock  wave.  This  agrees  with  the  observation  made  by 
Bazhenova  and  Soloukhin  (1959)  that  the  merging  of  the  flame  front  with  the 
incident  shock  wave  is  not  a  necessary  condition  for  the  detonation  wave 
formation.  The  simulations  also  show  a  large  pressure  overshoot  when  the 
forward  moving  wave  coalesces  with  the  incident  shock  wave.  This  has  been 
observed  earlier  in  the  experiments  of  Urtiew  and  Oppenheim  (1966). 

The  numerical  simulations  presented  in  this  paper  show  that  the  one¬ 
dimensional  reactive  flow  model  with  detailed  chemical  kinetics  can  be  used 
to  elucidate  some  of  the  details  of  weak  ignition  behind  incident  shocks  and 
the  subsequent  transition  to  detonation.  The  predictions  of  the  model  are  in 
qualitative  agreement  with  experimental  observations.  However,  the  model 
does  not  include  multi -dimensional  phenomena  such  as  turbulence  and  boundary 
layer  growth  which  play  an  important  part  in  any  quantitative  study  of  the 
transiton  to  detonation.  Currently  a  two-dimensional  reactive  shock  model 
exists  but  it  uses  a  parameterized  model  for  energy  release  (Oran  et  al. , 
1981b;  Oran  et  al. ,  1982).  The  model  is  now  being  extended  to  include  a 
detailed  chemical  kinetic  scheme.  Calculations  with  this  new  model  would 
show  the  effects  of  transverse  waves  and  boundary  layers  on  the  transition  to 
detonation. 
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Thble  I.  H2-02  Elementary  Reaction  Mechanism 


Reaction 

ki  * 
A(b) 

AT®  exp  (-C/T)(a) 

B 

c(b) 

References^ 

H  +  HO  t  0  +  H, 

i.uo(-iu) 

1.00 

3.50(+03) 

til 

3.00(-l4) 

1.00 

4.48(+03) 

[i! 

H  +  HO-  t  H,  +  0, 

4.20(-ll) 

0.00 

3. 50( +02) 

Ill 

9.10(-11) 

0.00 

2.9l(+04) 

til 

H  +  H02  t  HO  +  HO 

4.20(-10) 

0.00 

9.50(  +02) 

til 

2.00(-ll) 

0.00 

2.02(+04) 

ID 

H  +  HO-  t  0  +  H,0 

8.30(-ll) 

0.00 

5.00(+02) 

(21 

1.75(-12) 

0.45 

2.84(+04) 

kr  «  kf/Kc 

H  +  H,0,  t  HO,  +  H, 

2.80( -12) 

0.00 

1.90(+03) 

Ill 

1.20{-12) 

0.00 

9.1*0(+03) 

Ill 

H  +  H,0,  t  HO  +  H,0 

5.28(-10) 

0.00 

U.50(+03) 

(U 

3.99(-10) 

0.00 

l*.05(+OU) 

kr  «  kf/Kc 

HO  +  H,  t  H  +  H,0 

1.83(-15) 

1.30 

1.8M+03) 

(31 

1.79(-lU) 

1.20 

9.6K+03) 

1 3 1 

HO  +  HO  t  H2  +  02 

1.09(-13) 

0.26 

1.4T(+04) 

kf  m  kr*c 

2.82(-ll) 

0.00 

2.42(+04) 

(41 

HO  +  HO  t  0  +  H,0 

1.00(-l6) 

1.30 

0.00(+00) 

[3 1 

3.20(-15) 

1.16 

8.77(+03) 

kr  «  kf/Kc 

HO  +  HO,  t  H,0  +  0, 

8.30(-ll) 

0.00 

5.03(+02) 

(51 

2.38(-10) 

0.17 

3.69(+OU) 

*r  "  VKc 

HO  ♦  H,0  $  HO,  +  H, 

1.70(-11) 

0.00 

9.10(  +02) 

(U 

U.70(-ll) 

0.00 

1.65(+04) 

(U 

HO  ♦  H,  J  HO  +  H,0 

1.20(-12) 

0.00 

9.1*l(+03) 

(41 

1.33(-1U) 

0.U3 

3.62( +04) 

*r  *  VKc 

H02  +  H0?  t  H?02  +  02 

3.00(-ll) 

0.00 

5.00(+02) 

[2! 

1.57(-09) 

-0.38 

2.20(+04) 

kr  *  VKc 

•Bible  I.  (Continued)  H2-02  Elementary  Reaction  Mechanism 


Reaction 

■  Al® 

A(b) 

exp  (-C/T)(a) 

B 

C<*> 

References1 

0+  HO  X  H  +  0, 

2.72(-12) 

0.28 

-8.10(+01) 

kf  *  krKC 

3.70(-10) 

0.00 

8.45(+03) 

fm  r  c 

0  +  H02  X  HO  ♦  02 

8.32(-ll) 

0.00 

5.03(+02) 

151 

2.20(-ll) 

0.18 

2.82(+04) 

kr  «  kf/K( 

o  +  h2o2  X  h2o  +  o2 

1. U0( —12) 

0.00 

2.12(+03) 

12! 

5.70(-lU) 

0.52 

4.48(+04) 

kr  *  VKc 

0  +  H,0,  t  HO  ♦  HO, 

1.4o(-12) 

0.00 

2.13(  +03) 

121 

2.07(-15) 

0.64 

8.23( +03) 

‘r  '  VK< 

H+H+MtH,  +M 

1.80(-30) 

-1.00 

0.00(+00) 

[11 

3.70(-10) 

0.00 

4.83(-04) 

HI 

H  +  HO  +  M  X  H,0  ♦  M 

6.20( -26) 

-2.00 

0.00(+00) 

HI 

5.80(-09) 

0.00 

5.29(+04) 

[11 

H  ♦  02  ♦  M  X  H02  ♦  M 

4.l4( -33) 

0.00 

-5.00(+02) 

111 

3.50(-09) 

0.00 

2.30(+0U) 

HI 

HO  ♦  HO  ♦  M  X  H,0,  ♦  M 

2.50(-33) 

0.00 

-2.55(+03) 

Ill 

2.00(-07) 

0.00 

2.29(+04) 

fll 

0  +  H  +  M  t  HO  ♦  H 

8.28(-29) 

-1.00 

0.00(+00) 

161 

2.33( -10) 

0.21 

5.10( +04) 

kr  -  kf/K( 

0  +  HO  ♦  M  X  HO,  +  M 

2.80(-3l) 

0.00 

0.00(+00) 

[61 

l.lO(-OU) 

-0.43 

3.22(  +04) 

kr  -  kf/K( 

O+O  +  MtO,  +  M 

5-20( —35) 

0.00 

-9.00( +02) 

HI 

3.00(-06) 

-1.00 

5.91»(+04) 

Ill 

(a)  Bimolecular  reaction  rate  constants  are  given  in  units  of  cm3/(aolecule  sec). 
7e  molecular  reaction  rate  constants  are  given  in  units  of  cm6/ (molecule 2  sec) 

(b)  Exponentials  to  the  base  10  are  given  in  parenthesis}  i.e. ,  1.00(-10)  * 

1.00  x  10*1°. 

(c)  The  references  are:  (l)  Baulch  et  al. ,  1972;  (2)  Haapaon  and  Garvin,  19T5; 

(3)  Oohen  and  Westberg,  1979.  (^)  Olson  and  Gardiner,  1977}  (5)  Lloyd,  1974} 
(6)  Bahn,  1968. 


Ikble  II.  fhrameters  for  the  Weak  Ignition  Study 


Ihraaeter 


Eriver  Section  Driven  Section 


length 

Gas  Mixture 

Initial  Temperature 

Initial  Pressure 

Incident  Stock  Velocity 

Tesperature  behind  Incident  Stock 

Pressure  behind  Ihcldent  Stock 

Chapaan-Jouguet  Detonation  Velocity 


30  cat  370  cat 

He  H2:02:U2/2:1sU 

298  K  298  K 

9  ata  0.1  atm 

l.l*xl05  ca/s 
918  K 
1.39  ata 
1.9X10S  ca/s 


Thble  III.  Time  History  of  the  Forward  Moving  Reaction  Wave. 


Has  Velocity  of  Wave  Fluid  Velocity  Sound  Velocity 

(us)  (ca/s)  (cm/s)  (cm/s) 


Ahead 

Behind 

Ahead 

Behind 

609.03 

2.T0(-K>5) 

1.09(405) 

1.30(+05) 

7.15(-K)U) 

1.10(405) 

618.63 

2.T2('K>5) 

1.08(405) 

1.50(405) 

7.15(+0l*) 

1.12(405) 

626.03 

2.75(405) 

1.08(405) 

1.60(405) 

7.10(>OU) 

1.15(405) 

631.83 

2.80(+05) 

1.08(405) 

1.85(405) 

7.08(+0k) 

1.1  U( 405) 

663.83 

3.20(+05) 

1.07(405) 

2.10(405) 

7.20(+0l») 

1.18(405) 

688.13 

2.60(+05) 

0 

2.10(405) 

U.OU(+OU) 

1.28(405) 

Note:  Exponentials  to  the  base 

10  are  given  In  parenthesis, 

,  l.e. ,  2.70(405)  -  2.70x10 

Thble  IV.  Time  History  of  the  Backward  Moving  Reaction  Wave. 


Time 

(us) 

Velocity  of  Wave 
(ca/s) 

Fluid  Velocity 
(ca/s) 

Ahead  Behind 

Sound  Velocity 
(ea/s) 

Ahead  Behind 

609.03 

-7.2(40U) 

1.08(405) 

8.3(405) 

7.2(401*) 

1.06(405) 

618.63 

-7.3(401*) 

1.08(405) 

7.0(40U) 

7.2(401*) 

1.11(405) 

626.03 

-7.5(40U) 

1.09(405) 

6.0(401*) 

7.2(401.) 

1.13(405) 

631.83 

-7.6(401*) 

1.09(405) 

5. 5(401*) 

7.2(401*) 

I.IM405) 

numerical  simulations 


TIME  Ip*) 


i 
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Figure  2.  A  position-time  diagram  of  the  main  events  occurring  in  the  shock 


tube  sinulation:  S,  incident  shock;  CS,  contact  surface;  RC, 
reactive  centers,  HW,  reaction  wave;  (a)  pre-ignition  regime,  (b) 
quasi-eteady  shock -react ion  conplex,  (c)  formation  of  reactive 
centers,  (d)  hot  spot  formation  leading  to  overdriven  detonation 
and  (e)  detonation  relaxation. 


time  =  135 


Figure  3.  The  spatial  variation  of  the  temperature  and  the  OH  mole  fraction 


TEMPERATURE  <K) 


4 


Figure  8.  The  pressure,  velocity  and  tenperature  distributions  across  the 
system  at  (a)  599*33  Ms  and  (b)  609*03  Ms. 
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